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1  chose  this  thesis  topic  desllng  with  reconflgurable  control 
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David  Potts  for  the  AFWAL/Flight  Dynamics  Laboratory  (FDL) . 
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and  then  used  output  feedback  to  design  a  control  law  for  the  A-7D. 

The  technique  used  to  design  the  control  law  was  developed  by  Professor 
Brian  Porter  at  the  University  of  Salford,  England. 
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This  theala  uaea  the  deaign  procedurea  developed  by  Profeaaor 
Porter  at  the  Unlveraity  of  Salford,  England  In  an  attempt  to  deaign 
a  multivariable  tracker  control  law  for  the  A-7D  Dlgltac  II  Aircraft. 

Some  of  the  llmltatlona  and  problems  associated  with  this  design  pro¬ 
cedure  are  uncovered  In  this  study. 

A  six-degree -of- freedom  aircraft  model  is  developed  and  Is  then 
modified  to  a  form  that  is  required  by  the  design  procedure.  The  theory 
used  for  the  design  determines  the  necessary  arrangement  of  the  equations. 
A  tracker  control  law  is  first  designed  for  one  flight  condition.  Then 
it  is  checked  for  robustness  by  applying' the  control  law  at  a  different 
flight  condition  and  also  by  removing  the  rudder  from  the  inputs.  A 
design  computer  program  called  MULTI  is  developed  to  perform  the  compu¬ 
tations  and  simulations. 

It  is  found  that  the  design  techniques  developed  by  Professor 
Porter  are  valid,  but  that  they  are  not  applicable  to  all  systems.  A 
problem  occurs  when  the  inputs,  as  with  an  aircraft,  are  bounded. 

Problems  may  also  be  encountered  when  the  sensor  and  actuator  models 
are  incorporated  into  the  design.  Therefore,  for  this  study,  the  sensor 
and  actuator  models  are  removed  and  approximated  as  unity. 

More  work  is  needed  in  this  research  to  expand  knowledge  about  the 
selection  of  the  adjustable  parameters  in  order  to  develop  a  better 
design  and  to  more  effectively  utilize  the  basic  design.  Further  work 


la  also  needed  to  prefect  the  useabillty  of  the  progran  MULTI.  This 
thesis  provides  s  good  stepping  atone  to  a  better  understanding  of 
this  design  technique  and  Its  applicability. 
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Design  of  e 

Multivariable  Tracker  Control  Law 
for  the  A-7D  Digitac  IX  Aircraft 


I.  Introduction 

As  technology  continues  to  advance,  the  use  of  digital  flight 
control  systems  on  aircraft  grcws.  Future  aircraft  will  have  the 
normal  set  of  primary  control  surfaces  (ailerons,  spoilers,  flaps,  and 
horizontal  stabilizers)  split  into  Independently  controlled  sections. 

This  thesis  attempts  to  develop  a  tracker  control  law  with  the  primary 
surfaces  split  in  this  fashion.  In  order  to  learn  the  techniques  for 
implementing  the  theory,  the  decision  was  made  to  recombine  some  of  the 
surfaces  in  order  to  arrive  at  a  model  that  was  more  manageable.  The 
resulting  model  has  only  one  split  surface  (the  horizontal  stabilizer). 

Background 

IE  had  been  recently  suggested  chat  It  is  possible  to  design  a  direct 
digital  flight  control  system  that  offers  noolnteractlng  control  of  various 
flight  modes  which  allows  the  aircraft  Co  maneuver  in  ways  that  are  not 
possible  using  conventional  techniques.  Bradshaw  and  Porter  (Kef.  3) 
have  developed  a  synthesis  method  for  this  type  of  control  which  uses 
fast-sampling  error -actua tad  digital  control.  The  implementing  equations 
necessary  to  perform  this  design  have  been  developed  by  Porter  under 


contract  to  the  Air  Force  (Ref,  12).  This  thee  la  took  these  equations 
end  this  theory  and  used  them  to  investigate  a  flight  control  system 
for  the  A-7D  aircraft. 

Problem 


The  object  of  this  thesis  Is  to  develop  a  comprehensive  aircraft 
model  for  the  A-7D  aircraft  and  then  to  use  Professor  Porter's  theory 
to  design  a  tracker  control  law.  The  result  Is  a  review  of  this  new 
design  technique,  and  a  commentary  on  its  limitations,  problem  areas, 
and  performance. 


A  tracker  control  law  is  designed  in  order  to  gain  a  complete  under¬ 
standing  of  Porter's  methods.  The  resulting  design  is  Chen  checked  for 
robustness  in  two  ways: 

1.  The  control  law  Is  used  at  a  different  flight  condition 
to  see  how  it  responds  over  a  vide  range  of  parameter 
changes . 

2.  The  rudder  Is  made  inoperative  and  the  control  lav  used 
to  see  how  this  aircraft  disability  is  handled. 

Approach 

This  thesis  is  limited  to  an  investigation  of  the  theory  developed 
by  Professor  Porter  (Ref.  4).  The  resulting  design  is  not  claimed  to  be 
the  best  possible,  but  rather  the  best  the  this  engineer  could  design  In 
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the  tine  allowed.  The  design  la  cheeked  for  robustness,  and  the  results 
are  discussed.  The  entire  design  process  Is  done  using  an  Interactive 
computer  package  called  MULTI.  MULTI  has  been  developed  by  Captain  Doug 
Porter,  Lieutenant  Joseph  Smyth,  and  this  author.  The  progran  utilizes 
the  equations  developed  by  Professor  Porter.  The  steps  taken  In  completing 
this  thesis  are: 

1.  Develop  the  necessary  aircraft  equations  at  three  flight 
conditions. 

2.  Develop  the  interactive  computer  software  package  MULTI 
that  performs  the  design  and  simulations. 

3.  Design  a  tracker  control  law  using  the  appropriate  method. 

4.  Check  the  resulting  design  for  robustness. 

Each  step,  by  itself,  could  have  been  expanded  Into  a  separate 
thesis.  Since  there  are  three  possible  approaches  to  the  problem  (see 
Chapter  III),  step  3  proved  to  be  very  difficult.  Numerous  trials  are 
attempted  before  an  acceptable  design  that  works  Is  finally  reached. 

Assumptions 

The  assumptions  made  In  this  thesis  deal  with  the  aircraft  model, 
and  the  aircraft  operation.  First,  the  aircraft  Is  assumed  to  make  only 
small  perturbations  about  a  trimmed  straight  and  level  flight  condition. 
Secondly,  the  aircraft  equations  of  motion  developad  are  In  the  stability 
axis  reference  frame.  Thirdly,  the  coefficients  of  the  control  derivatives 


are  normally  for  a  pair  of  confcrol  surfaces,  are  halved,  given  the 
proper  sign,  and  used  whenever  the  control  surface  Is  spilt  Into  two 
Independent  sections. 

Presentation 

This  thesis  Is  composed  of  7  chapters.  Chapter  II  develops  the 
aircraft  models  and  Appendix  A  provides  the  details  and  the  equations 
used.  Chapter  III  presents  the  theory  used  for  the  design  of  the  digital 
control  law.  Three  approaches  are  discussed,  as  well  as  when  to  use 
each  one.  Chapter  IV  discusses  the  development  of  the  computer  program 
MULTI  which  is  needed  to  compute  the  design  and  simulate  the  results. 

The  program  handles  all  three  design  approaches.  Chapter  V  describes 
several  approaches  taken  in  an  attempt  to  find  an  approach  that  produces 
valid  results.  Uhen  more  than  one  approach  can  be  used,  it  is  necessary 
to  find  out  which  one  Is  better  and  then  to  use  that  one.  Chapter  VI 
discusses  the  results  obtained.  The  results  are  collected  in  Appendix  B 
for  easy  location  and  viewing.  Chapter  VII  presents  the  conclusions  and 
recommendations . 
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II.  The  Aircraft  Models 


Introduction 

The  model  for  the  A-7D  aircraft  presented  by  Potts  (Ref.  15)  is 
revised  and  a  new  model  is  developed.  Using  the  same  six  degree-of- 
freedom  (D-O-F)  aircraft  model  and  the  same  approach  as  presented  by 
Potts  (Ref.  15).  a  "modified "model  is  derived.  The  model  considers 
many  control  surfaces  to  be  Independent  which  are  traditionally  con¬ 
sidered  to  act  as  one  unit.  For  example,  the  horizontal  stabilizer  is 
divided  into  two  Independent  control  surfaces,  the  right  and  left  hori¬ 
zontal  stabilizers.  A  linear  model  is  obtained  by  linearizing  the  non- 
flnear  six  degree-of-freedom  equations  about  a  normal  operating  point. 
After  a  revised  model  Is  generated  for  the  first  flight  condition  at  0.6 
Mach  at  an  altitude  of  15,000  feet,  two  additional  flight  conditions  arc 
considered.  Models  are  also  developed  for  Mach  numbers  of  0.18  and  0.S 
altitudes  of  2,000  feet  and  35,000  feet  respectively.  The  model  for  0.3 
Mach  is  assumed  to  be  a  landing  configuration.  In  this  manner,  the  new 
control  laws  are  checked  for  robustness  over  a  wide  range  of  flight 
conditions. 

Aircraft  Description 


The  A-7D  is  a  single  seat,  light  attack  aircraft  with  moderately 


swept  wing  and  tall  surfaces,  and  la  powered  by  a  single  TF41-A-1 
engine  (Ref. 17).  The  describing  data  for  the  aircraft  as  found  In 
Reference  9  Is  given  in  Table  1. 


TABLE  I 

A-7D  Descriptive  Data 


Item 

Dimension 

Unit 

Fuselage 

Length 

45.4 

ft 

Wing  Area 

375 

ft2 

Wing  Mean 
Geometric  Chord 

10.8 

ft 

Horizontal 
Stabilizer  Area 

56.2 

ft* 

-Horizontal 

Slab  m.g.c. 

6.1 

ft 

Distance  from 

0.25  Wing  m.g.c. 
to  0.25  Hori¬ 
zontal  Slab 
m.g.c. 

16.2 

ft 

Weight 

25,238 

lbs 

Svtem  Models 

The  models  used  to  describe  the  A-7D  at  various  flight  conditions 
are  developed  using  the  same  6  D-O-F  aircraft  equations  given  by  Potts 
in  Reference  IS.  By  using  the  same  technique  as  Potts ,  control  of 
lateral  motion  is  made  possible  by  using  a  longitudinal  control  surface, 
and  control  of  longitudinal  motion  is  made  possible  by  using  a  lateral 
control  surface.  This  Is  essential  for  reconfiguration  to  be  possible 
after  the  loss  of  a  primary  control  surface.  The  details  of  this  technique 
are  found  in  Appendix  A.  The  control  surfaces  that  are  considered  as  in¬ 
put  are  the  ailerons,  the  spoilers,  the  flaps,  the  rudder,  and  the  hori¬ 
zontal  stabilizer.  The  ailerons,  spoilers,  and  the  horizontal  stabilizer 
are  divided  into  separate  individual  control  surfaces  (right  and  left), 
while  the  rudder  and  flaps  are  left  as  single  inputs.  In  this  manner  there 
are  a  total  of  8  inputs  considered.  There  are  likewise  8  outputs.  These 
include  the  forward  velocity  u,  the  flight  path  angle  y ,  the  pitch  rate  q, 
the  pitch  angle  6,  the  roll  rate  p,  the  roll  angle  0,  and  the  yaw  rate  r. 
The  system  is  square,  that  is,  the  number  of  inputs  is  equal  to  the  number 
of  outputs.  By  including  the  individual  control  surfaces,  the  equations  of 
motion  are  not  decoupled  to  separately  describe  lateral  and  longitudinal 
motion. 

The  stability  derivatives  that  are  needed  in  the  equations  of  motion 
are  given  in  Appendix  A.  The  derivatives  that  are  not  given  in  either 
Reference  1  or  Reference  9  are  computed  by  conventional  aerodynamic 
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techniques  with  the  eid  of  e  computer  package  called  digital  DfcTOOM 
(Ref.  10). 


Sumnary 

This  chapter  gives  a  description  of  the  A-7D  aircraft  and  identifies 
the  control  surfaces  which  are  considered  as  Inputs  for  the  6  D-O-F 
models  developed.  Models  for  three  flight  conditions  are  developed  with 
0.6,  0.18,  and  0.8  Mach  at  altitudes  of  15,000,  2,000  and  35,000  feet, 
respectively.  With  the  models  developed,  the  next  step  Is  to  present  the 
multivariable  control  theory  by  which  the  new  control  laws  are  designed. 


III.  Multivariable  Digital  Control  Theory 


Introduction 

The  design  techniques  in  this  thesis  for  development  of  a  tracker 
control  lav  use  the  theory  developed  by  Professor  Brian  Porter  (Ref.  3). 

The  concepts  Involved  and  the  describing  equations  are  presented  in  this 
chapter.  Only  the  digital  approaches  are  discussed,  although  analog 
control  laws  based  on  the  same  concept  are  possible.  These  techniques  are 
used  to  synthesize  a  samp led -data  hybrid  control  system  which  consists  of 
a  continuous- time  plant  and  a  digital  controller  that  produces  a  control 
input  signal  which  is  piece-wise  constant 'for  each  sampling  period. 

Professor  Porter  has  developed  three  design  procedures.  All  three 
of  the  procedures  are  presented  in  this  chapter.  The  procedure  to  use 
is  determined  by  the  first  markov  parameter.  The  first  markov  parameter 
is  defined  as  the  product  CB  and  is  discussed  in  this  chapter. 

Unknown  Plant  Design 

In  many  cases,  a  state  equation  model  is  not  available  to  describe 
a  system.  If  this  is  the  case,  the  steady-state  transfer  function  matrix 
G(0)  can  be  determined  from  off-line  tests  performed  on  the  system  (Ref.  13) 
provided  that  the  plant  is  stable.  If  the  state  equations  are  available, 
then  G(0)  can  be  determined  using  Equation  (1)  for  the  system  represented 
by  Che  continuous -time  state  and  output  equations,  &  *  Ax  +  Bu  +  Dd 

and  y  -  Cx  +  Du. 
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0(0)  -  -C  A'1  B  (1) 

where  C(0)  £  R*x  m 

In  either  case  G(0)is  the  transfer  function  matrix  G(X) ,  as  defined  in 
Equation  (2),  with  X"  0,  that  is, 

GOO  -  C  UIn  -  A  )_1  B  (2) 

where  n  ■  number  of  states 

In  order  for  integral  action  to  preserve  stability,  G(0)  must  have 
rank  equal  to  the  number  of  outputs.  This  condition  also  requires  that 
the  the  number  of  outputs  be  less  than  or  equal  to  the  number  of  Inputs 
(  X  £  si).  This  procedure  assumes  that  the  eigenvalues  of  the  open 
loop  plant  matrix  A  €  r"310  lie  in  the  open  left  half-plane  (Ref .  12  ). 

Digital  error-actuated  controllers  are  described  by  state  and  output 
equations  of  the  form  (Ref.  13  ) 

V,  ■  +  H 

*k  ■  r*k 

where 

xk  -  x(kT)  £  Rn 

uk  ■  u(kT)  £  Rm 

yk  -  y(kT)  e  R* 

-  d(kt)  £  RP 

-Tk  -  «xp  (AT) 

-  CT  exp  (At) 


Ld 


(3) 

(4) 


(state  vector) 

(control  input  vector) 

(output  vector) 
(disturbance  vector) 


B  dt 
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A  ■  J  np  (At)  D  dt 
uo  ” 

r-  c 

In  these  equations  T  is  the  sampling  period. 

When  error>actuated  digital  control  Is  applied  to  plants  described 
by  Equations  (3)  and  (4),  the  proportional  plus  Integral  control  law 
equation  that  results  is: 


uk  “  °<rTt£*ek  +  £™**k  (5) 

where 

efc  •  e(KT)  •  6.  (error  vector) 

l 

vk  ”  ^k(KI)  G  R  (command  input  vector) 

_  a 

-  a^KT)  €  R  (Integral  of  error) 

*k+l  "  *k  +  Tek 


In  Equation  (5),0<ris  a  design  parameter  that  determines  the  desired  ratio 
of  proportional  to  Integral  action.  £  is  a  normalizing  factor,  and  vk  is 
the  command  input  vector.  If  ©<r,  T,  and  K  are  chosen  so  that  all  the 
closed-loop  eigenvalues  lie  within  the  digital  unit  disc,  then  jLim  e  -  0, 
and  set-point  tracking  occurs  (Ref.  13).  According  to  reference  13  the  para¬ 
meter  K  is  determined  by  the  relationship 

K  -  GT(GG  V1  Z.  (6) 


where 


In  Equation  (6),  X  la  a  matrix  that  la  used  aa  an  additional  "tuning" 
parameter  for  better  reaulta.  Figure  1  la  a  block  dlagraai  for  thla 
control  scheme. 


Regular  Planta 

Planta  that  are  described  by  state  and  output  equations 


i  -  Ax  +  Bu 
y  »  Cx 


(7) 

(8) 


can  be  transformed  by  use  of  a  transformation  matrix*  if  necessary,  so 
that  the  new  state  and  output  equations  have  the  form 
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where 


€ 

€. 

£ 

<= 

£ 

£ 

£ 

£ 

£ 

£ 

£ 


R(n-/)  x  (n-f) 
R(n-f) 

/*  <«»-*> 

4x4 

R 

lx  l 

R 

R^*  (n-4) 

lx/ 

R 


The  first  Markov  parameter  is  defined  as  *or  equations  of  this 

form.  If  the  rank  of  £2-2  "  ^  '  than  the  plant  is  referred  to  as 

"regular".  For  plants  of  this  type,  the  governing  control  law  equation 
Is  (Ref.  3  ). 

O(kT)  fil^RT)  +  KjCZ(kT)  3  <U> 


where 


So* 


lx  l 


R 
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Por  this  control  lav,  where  e(kT)  »  v(kT)  -  y(kT) f  the  output  vector 
y(kT)  tracks  any  constant  command  vector  v(kT)  (Ref.  3).  The  relation¬ 
ship  between  and  Is 

-  <*.%  (12) 

where  o^.is  an  adjustable  design  parameter.  Is  defined  by  the  equation 

"  <C2B2)'l*Z  (13) 

where 

%  o-j . 

-1  <  <  1  -  Oj  )  <  1 
03  €.  R 

The  closed-loop  system  responds  faster  and  has  less  Interaction  as  T 
Is  decreased  (Ref.  3).  For  systems  of  this  form.  Equation  (9), the 
transmission  zeros  are  the  set  derived  from  the  relationship 

*T  ■  (|  -n  -  !„  -  “u  +  “uc2lcl  I  ’  °J  «*> 

The  transmission  zeros  must  lie  within  the  unit  disc  for  complete  stability 
as  T  -v  0.  Transmission  zeros  that  lie  on  the  unit  disc,  or  outside  of  It, 
cen  ceuse  instebilitles.  Also,  if  paths  to  transmission  zeros  cross  Into 
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unittbU  regions,  Cher*  exists  a  range  of  T  choc  produces 
Instability.  For  tracking  to  occur,  all  of  the  closed- loop  poles  must 
lie  within  the  unit  disc.  Figure  1  Is  also  s  block  diagram  for  this 
type  of  control  scheme . 


Irregular  Plants 


For  high-performance  systems  described  by  Equation  (9)  and  Equation 
(10),  if  the  first  Markov  parameter  has  rank  less  than  ,  then  the 
plant  is  said  to  be  “Irregular".  It  is  assumed  that  Is  full  rank  in 
such  systems  so  that  the  rank  deficiency  Is  within  the  matrix  £2.  If  this 
is  the  case.  It  Is  necessary  to  Introduce  a  vector  of  extra  Masureamnts 
U(t)  (Ref.  4  )  such  that 

h«>  -  Or  1.  £ J 


where 


Fj  €  1 


Ax  (n-jf) 


u  £  * 


Ax  A 


The  addition  of  the  extra  measurements  changes  the  form  of  e(t)  so  that 
e(t)  •  v(t)  -  W(t).  Fj  and  F^  are  defined  by  the  equations 


Si  ♦  s*n 


£2  ♦  ma12 


(16) 

(17) 


Za  Equations  (16)  and  (17),  M  is  the  transducer  amtrlx  of  dlsmnslon 


for  constant  command  input a,  ao  that  in  ateady  atate  H(t)  -  y(t)  (Ref.  4). 
Therefore,  in  the  ateady  atate 


Jim^  e(kT)  -  lio^  v(kT)  -  W(kT)  -  0  (19) 

This  ensures  tracking  of  the  cosaaand  vector  in  steady  state  operation. 

Proper  choice  of  M  "makes"  the  plant  regular.  This  requires 
that  FjBj  have  rank  X.  .  Since  is  assumed  to  have  rank  X  at  the  start, 

M  is  chosen  so  that  (Ref.  4  ). 

rank  F2  -  rank  (  C2  +  M  Al2  )  -  X  (20) 

Once  the  plant  is  "made”  regular  by  proper  choice  of  M  ,  the 
governing  control  law  is  equation  (11)  where  now 

5o  "  (  -2-2  >’1*  I  (21) 

The  transmission  zeros  are  the  set  derived  from  the  relationship 


*T  "  [  V/  -  T  AU  +  TA12^21fi(  "  °]  <22> 


The  seme  ideas  that  applied  to  regular  plants  apply  to  irregular  plants 
once  M  is  chosen  so  as  to  "stake"  the  irregular  plant  appear  regular. 
Figure  2  is  a  block  diagram  for  this  control  scheme. 
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Irregular  Plant  Control  Diagram 


Suwaary 


This  chapter  presents  the  three  design  concepts  developed  by  Professor 
Porter.  For  further  details,  the  reader  is  encouraged  to  see  the  references 
given  in  this  chapter.  The  three  categories  of  plants  are:  unknown, 
regular  (rank  CjBj  *  X  ) »  *nd  irregular  (rank  <  X  )  •  The  next  chapter 
presents  the  coeiputer  package  developed  to  employ  these  concepts.  Chapter 
1  shows  how  the  proper  design  procedure  is  chosen  and  applied. 


Introduction 


MULTI,  an  interactive,  user-oriented  computer  program,  is  devel¬ 
oped  to  design  and  simulate  the  control  laws  for  unknown,  regular,  and 
irregular  plants.  The  computer  program  is  written  in  FORTRAN  V  code. 

It  is  developed  along  the  lines  of  a  computer  package  called  PAK200 

f 

which  was  received  from  the  University  of  Salford,  Department  of  Aero¬ 
nautical  and  Mechanical  Engineering  (Ref.  6  ).  The  program  utilizes 
digital  theory  as  opposed  to  continuous  methods.  The  MULTI  computer 
package  contains  approximately  40  options  which  give  the  user  an  inter¬ 
active,  iterative  approach  in  the  design  and  simulation  of  control  laws 
for  linear,  multivariable  plants.  The  control  law  assures  that  the  out¬ 
put  tracks  the  input  and  that  disturbance  rejection  is  accomplished. 

This  chapter  relates  the  requirements  that  are  set  for  the  pro¬ 
gram,  the  program  design  emphasis,  the  constraints  and  problems 
faced  in  the  development  of  MULTI,  and  the  actual  program  structure 
of  MULTI.  A  User's  Manual  for  MULTI  and  a  Programmer's  Manual  for  MULTI 
are  found  in  Captain  Doug  Porter's  thesis  (Ref.  14  )  along  with  a  complete 
listing  of  the  program. 

Program  Requirements 

The  requirements  for  the  computer  program  are  set  so  that  MULTI 
can  take  full  advantage  of  the  techniques  used  in  designing  the 
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discrete-time  tracking  systems  for  linear  multivariable  plants. 
These  requirements  include: 


li  The  computer  package  must  be  fully  interactive,  user- 
oriented,  and  allow  for  an  interactive  design  process. 

2.  The  program  must  allow  the  user  to  input  data  from  the 
terminal  keyboard  or  a  data  file  and  store  pertinent 
data  in  local  files  upon  normal  program  termination. 

3.  The  package  should  include  design  capability  for  unknown, 
regular  and  irregular  plants. 

4.  The  package  must  be  able  to  recognize  when  singular 
matrices  are  formed  and  direct  the  user  to  apply  an 
alternate  design  technique. 

5.  The  package  must  be  able  to  form  a  measurement  matrix  from 
terminal  input. 

6.  The  package  must  include  a  discrete -time  simulation  in 
order  to  evaluate  the  control  laws  developed. 

7.  The  package  must  include  a  plotting  capability. 


These  requirements  result  in  a  software  package  that  is  very  versatile. 


Program  Design  Emphasis 

When  a  computer  program  is  divided  into  several  modules,  the  content, 
of  each  module  and  the  interconnections  between  the  modules  can  signi¬ 
ficantly  affect  the  operation  and  complexity  of  the  resulting  program. 
Coupling  and  cohesion  are  two  qualities  that  are  used  to  check  the  over¬ 
all  design. 

Coupling  is  the  measure  of  the  strength  of  the  interconnections 
between  one  module  and  another.  It  is  essential  that  coupling  be  kept 
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at  a  minimum  so  that  changes  In  one  module  do  not  severely  affect  the 
other  modules.  When  coupling  is  retained  at  a  low  level,  the  ease  of 
finding  and  correcting  program  bugs  is  enhanced  and  the  program's  life 
is  increased  since  quick  modifications  are  made  possible.  Labeled 
comnon  blocks  are  used  rather  than  blank  common  blocks  so  that  the  level 
of  coupling  is  reduced. 

Cohesion  is  the  degree  of  functional  relatedness  of  processing 
elements  within  a  single  module.  It  also  has  a  direct  effect  on  the 
ability  of  a  programmer  to  maintain  and  modify  a  program.  A  high  level 
of  cohesion  is  desirable  and  is  kept  by  putting  only  functionally  re¬ 
lated  processes  within  each  module. 

In  MULTI,  the  major  functions  of  the  program  are  directed  to  lower- 
level  modules.  Every  module  is  designed  to  use  a  minimum  of  memory  core, 
and  every  labeled  common  block  is  designed  to  hold  only  those  elements 
needed  for  a  specific  purpose.  Through  this  type  of  design,  MULTI  is 
able  to  use  minimum  computer  memory  core,  have  minimum  coupling  and  retain 
a  high-level  of  cohesion. 

Design  Constraints 

There  are  two  main  design  constraints  when  developing  a  computer 
program  like  MULTI.  One,  the  computer  language  that  is  used,  and 
two,  the  availability  of  computer  core  memory  that  i6  allowed  when  using 
an  INTERCOM  system. 


FORTRAN  V  is  chosen  as  the  language  for  MULTI  since  it  is  the 


latest  version  of  FORTRAN  eve liable  on  the  ASD  CYBER  computer  system. 

Using  FORTRAN  V  does  not  restrict  the  use  of  the  ASD  Library  sub¬ 
routines,  even  though  they  are  coded  In  FORTRAN  IV,  since  all  are 
readily  compiled  into  the  more  powerful  version. 

Current  Interactive  computer  users  ere  constrelned  to  operete  in 
6SK  words  of  core  memory,  or  less.  Therefore,  MULTI  Is  designed  to 
operate  within  this  restriction.  The  core  memory  restriction  Is  met 
through  the  modular  design  of  the  program. 

Developmental  Problems 

The  original  program  design  is  edopted  from  e  computer  package 
called  PAK200  which  was  written  for  controller  design  simulation  at 
the  University  of  Salford,  England.  It  was  originally  conveyed  that 
the  development  of  a  similar  computer  peckage,  compatible  with  the  ASD 
CYBER  computer,  would  merely  Involve  the  substitution  of  certain  computa¬ 
tional  subroutines  from  the  ASD  Subroutines  Into  PAK200.  The  MULTI 
computer  package  is  now  a  completely  Interactive  design  and  simulation 
tool  for  all  three  types  of  plents.  The  main  body  of  the  code  from  PAK 
200  Is  entirely  contained  in  one  option  of  MULTI.  TWo  subroutines,  dealing 
with  the  differential  equation  formation  and  the  value  of  the  output  at 
the  end  of  each  time  increment,  are  also  retained  to  provide  compatibility 
with  the  original  simualtlon  code. 

MULTI's  development  faced  several  other  problems.  The  use  of  the 
FORTRAN  V  lenguage  presented  e  barrier  which  was  not  eliminated  until 
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several  months  Into  MULTI '•  design.  This  problem  Involved  the  use  of 
ehsrseter  statements,  first  available  in  FORTRAN  V,  which  produced 
sporadic  errors  during  MULTI's  many  modifications.  This  problem  is 
known  to  FORTRAN  V  subscribers  but  not  to  it*  users. 

The  constant  breakdown  of  computer  operations  is  a  big  problem. 

Just  when  it  seems  progress  is  being  made,  the  ASD  CYBER  system  would 
go  down.  This  caused  the  proposed  time  for  program  development  to  be 
increased  by  as  much  as  a  month.  This  is  a  problem  that  should  be 
considered  by  anyone  developing  or  modifying  a  program. 

The  final  version  of  MULTI  is  the  best  that  can  be  obtained  by 
control  engineer  programmers  within  the  time  constraints  imposed  on  the 
development,  testing,  and  simulation  of  a  model,  controller  and  computer 
package  for  the  newest  techniques  of  the  multivariable  design  of  this 
thesis.  MULTI  meets  all  the  specific  requirements  as  outlined  above. 

Program  Structure 

Since  MULTI  must  be  able  to  operate  in  65K  words  of  core  memory, 
its  program  structure  consists  of  a  group  of  various  modules.  Information 
between  these  modules  is  passed  by  using  the  data  base  concept  of  global 
storage.  To  allow  MULTI  to  be  fully  interactive,  the  program  uses  a  very 
simple  Program  Control  Interface. 

To  reduce  the  amount  of  memory  core  required  for  operation,  MULTI 
has  a  structure  consisting  of  several  modules,  called  overlays.  Each 
overlay  is  an  executable  program  and  combines  with  all  other  overlays  to 
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fora  a  .complete  functioning  overlay  structure.  There  is  one  main  over- 
lay  in  MULTI.  The  main  overlay  la  an  Executive  Program  directing  the 
actual  functions  of  the  program  to  13  lower-level  primary  overlays. 

This  overlay  structure  is  responsible  for  MULTI' s  operation  within 
the  defined  65K  memory  core  limit.  The  main  overlay  remains  in  executable 
memory  at  all  times.  As  each  computational  or  functional  requirement  of 
the  program  is  needed,  the  overlay  designed  to  perform  that  procedure  Is 
loaded  into  executable  memory  behind  the  main  overlay.  When  the  procedure 
Is  finished,  the  used  overlay  is  saved,  and  the  overlay  needed  for  the 
next  procedure  is  loaded  behind  the  main  overlay. 

Data  information  is  passed  between  overlays  in  MULTI  by  the  use  of 
labeled  comnon  blocks.  If  a  program  variable  in  an  overlay  is  not  listed 
in  a  labeled  comnon  block,  its  value  is  not  retained  when  the  next  over¬ 
lay  is  loaded. 

A  complete  description  of  overlays  and  how  they  are  used  in  MULTI  can 
be  found  in  the  MULTI  Programmer's  Manual  (Ref.  14  ) 

MULTI  uses  three  types  of  data  storage  methods  which  comprises  its 
data  base.  These  three  types  of  storage  methods  are:  local  storage, 
global  storage,  and  mass  storage.  Mass  storage  in  MULTI  entails  only  the 
use  of  sequential-access  files. 

Local  storage  la  the  storage  method  which  is  used,  for  program  variables 
that  are  needed  only  during  the  execution  of  a  particular  overlay.  As 
mentioned  above,  these  locally  stored  variable  values  are  lost  when  the 
overlay  is  finished  executing.  MULTI  has  several  matrices  which  are 
generated  during  program  execution  whose  values  are  destroyed  when  they 


are  no  longer  needed. 

Global  storage  is  used  for  program  variables  whose  element  values 
are  required  to  be  passed  between  overlays.  All  of  MULTI* s  global 
storage  is  accomplished  using  labeled  common  blocks.  The  variables 
defined  In  each  common  block  are  purposely  chosen  to  keep  each  over¬ 
lay  memory  core  requirement  at  a  minimum,  and  still  allow  for  the  necessary 
data  transfer  between  overlay  loading. 

Sequential-access  files  are  used  in  MULTI' s  external  data  input 
functions  and  in  MULTI's  mass  storage  capability.  A  sequential-access, 
data-lnput  file  can  be  built  by  any  user  for  use  with  the  program.  The 
specifics  and  proper  formatting  of  these  files  is  described  in  detail  in 
the  user's  manual  (Ref.  14).  MULTI  automatically  creates  two  separate 
memory  files  when  the  user  terminates  program  operation.  A  file  called 
"1CM0"  is  generated  to  hold  all  data  describing  the  plant  of  Interest, 
and  a  file  called  "MEM10"  is  created  to  hold  all  design  parameters. 

MULTI's  interface  design  is  similar  to  the  Interface  found  in  TOTAL 
(Ref.  g  ),  but  is  not  nearly  as  complex.  All  MULTI  operations  are  con¬ 
trolled  by  having  the  user  choose  various  option  numbers  which  correspond 
to  desired  computational  functions.  Since  there  is  minimal  input  data 
verification,  the  user  must  adhere  strictly  to  the  limitations  noted  in 
various  parts  of  the  MULTI  User's  Manual.  Every  completed  option  sets 
an  internal  "flag"  allowing  entry  into  other  options.  There  is  a  definite 
pattern  for  successful  program  operation  which  corresponds  to  a  normal 
design  process.  The  program  is  designed  to  require  the  user  to  flew  from 
data  input,  to  controller  design,  and  then  to  simulation.  Each  part  of 


the  design  cen  be  re -accomplished  Interactively,  but  the  user  la  not  per¬ 
mitted  to  proceed  in  the  flat*  until  each  preceedlng  design  step  is  accom¬ 
plished.  All  options  are  explained  in  the  MULTI  User's  Manual.  Data  access 
in  MULTI  is  accomplished  by  program  requests  and  by  the  use  of  a  special 
option  number  range.  The  user  can  view  certain  pertinent  matrix  combina¬ 
tions  and  other  data  of  interest  responding  affirmatively  when  the  program 
prints: 


ENTER  "0"  TO  OBTAIN  DATA  PRINTOUT 
ENTER  "1"  TO  SKIP  DATA  PRINTOUT... 

In  some  cases,  the  data  cannot  be  accessed  again,  since  it  is  contained 
on  a  local  storage  device,  unless  the  user  re-enters  the  same  option  with 
the  same  plant  data  and  design  specifications. 

The  second  method  of  data  access  is  by  the  use  of  options  having 
the  value  of  100  or  greater.  If  data  has  been  entered  using  OPTION  #1 
the  same  data  can  be  accessed  for  verification  by  using  OPTION  #101. 

This  method  of  data  access  is  developed  in  the  last  few  versions  of  the 
program.  The  original  in-line  format  of  displaying  data  takes  too  much 
time  during  the  complete  design  process.  The  final  data  access  method 
allows  for  easy,  selective  data  access. 

After  data  access  the  user  may  determine  that  values  are  not  entered 
correctly,  or  for  some  reason  need  to  be  changed.  The  user  can  use  one 
of  the  two  means  to  change  the  data.  The  option  related  to  the  data  in 
question  can  be  re-entered  or  the  user  can  terminate  MULTI  and  edit  the 
program- created  data  files.  This  last  means  of  changing  data  is  of 
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special  interest  when  changing  element  values  of  the  plant's  A,  B,  C, 
and  D  matrices. 

The  error  detection  and  recovery  aspects  of  MULTI  are  limited. 
Since  the  main  emphasis  of  the  thesis  was  not  originally  perceived  to 
entail  the  complete  design  of  a  computer  simulation  package,  MULTI's 
design  is  not  based  upon  the  ability  to  protect  the  user  from  input 
errors.  MULTI  is  designed  to  alert  a  user,  and  recover,  when  program 
flow  is  not  in  the  proper  order  or  when  certain  plant  deficiencies 
require  the  use  of  a  certain  controller  design. 

Input  errors  such  as  character  entries  for  numerical  data,  and 
vice  versa,  are  all  detected  by  the  CYBER  error  detection  capabilities. 
In  these  cases  the  user  is  allowed  to  re-enter  the  proper  information 
and  continue. 

MULTI  restricts  the  user  to  a  10  state,  10  input,  10  output  system 
with  2nd  order  actuators  and  sensors.  Simultaneous  simulations  can  be 
run  for  up  to  two  different  sampling  times  while  plots  can  be  generated 
which  display  up  to  four  superimposed  curves.  When  the  user's  requests, 
exceed  these  limitations,  the  entire  MULTI  program  is  aborted.  It 
is  suggested  that  any  MULTI  user  become  completely  familiar  with  the 
user's  guide  for  the  program. 

When  MULTI  is  aborted  the  CYBER  system  relates  the  nature  of 
the  error  which  causes  the  abnormal  termination  of  the  program.  It 
mist  be  stated  that  MULTI  is  not  as  dynamic  in  its  ability  to  receive 
input  data  as  other  programs  with  which  the  user  may  be  familiar  (l.e. 
TOTAL  (Ref.  8  ),  CESA  (Ref.  7  )).  At  present,  character  inputs  to 


display  current  matrix  values  or  system  design  values  are  not  available 
except  by  requesting  the  proper  option  number.  Also  at  present,  enter¬ 
ing  a  symbol,  rather  than  actual  data,  to  abort  an  option.  Is  not 
recognized  unless  specifically  noted  In  the  program  instructions  . 

Summary 

This  chapter  discusses  the  development  of  the  computer  program  MULTI. 
The  program  structure,  requirements,  and  design  is  presented  so  that  the 
reader  can  understand  what  is  done,  not  necessarily  how.  The  program  has 
areas  that  need  improvement,  but  these  can  only  be  identified  by  user 
involvement.  The  program  should  be  modified  before  release  so  that  it 
performs  as  desired  by  control  designers.  Further  details  of  the  program 
are  contained  in  Captain  Doug  Porter's  thesis  (Ref. 14  )  since  he  aided  in 
developing  the  package.  The  next  chapter  discusses  the  design  that  is 
performed  using  MULTI. 


V.  Design  Approaches 


Introduction 

This  ehspter  discusses  the  several  approaches  that  are  taken  to 
design  the  tracker  control  lav.  Since  there  are  three  design  procedures 
(see  Chapter  III),  the  first  thing  that  needs  to  be  determined  Is  the 
design  procedure  to  be  used.  Then,  once  the  proper  choice  Is  made,  the 
design  is  performed  using  the  program  MULTI. 

Looking  at  the  model  for  Mach  0.6  at  an  altitude  of  15,000  feet 
(see  Appendix  A)  It  is  obvious  that  the  B  matrix  does  not  have  full 
rank.  The  first  approach  attempted  is  the  use  of  the  psuedoinverse  of  B 
in  conjunction  with  the  regular  design  procedure  to  generate  a  control 
lav.  The  second  attempt  involves  recombining  some  of  the  Input  surfaces 
and  re-arranglng  the  equations  so  that  they  can  be  put  In  the  form  of 
Equation  (9)  and  Equation  (10)  In  Chapter  III.  Once  the  equations  are 
In  the  necessary  form,  the  unknown  design  procedure  Is  attempted  first. 
However,  §.(0)  proves  to  be  renk  deficient  so  the  regular  design  procedure 
Is  used  since  C?B?  had  rank  ■  X  »  The  first  choice  of  outputs  leads  to  a 
set  of  transmission  zeros  at  the  origin  and  an  uncontrollable  mode.  A 
second  choice  of  outputs  is  selected  with  Cjj^  rank  deficient.  Using  the 
Irregular  design  procedure,  a  measurement  matrix  M  is  selected  that  places 
the  transmission  zeros  at  -3.  This  gives  a  stable  response  and  a  control 
law  that  responds  to  a  desired  input  command  vector.  All  of  these  attempts 
are  explained  In  detail  in  this  chapter.  The  results  of  the  /Inal  attempt 
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arc  presented  to  show  the  resulting  tracker  control  law. 


Psuedolnverse  Approach 

The  theory  developed  by  Professor  Porter  considers  only  the  case 
for  n*>/t  (number  of  states  >  number  of  outputs).  However,  the  model 
presented  in  Appendix  A  has  n  ■  X  •  If  n  m  £,  then  the  B  matrix  Is 
equal  to  since  the  number  of  aero  rows  above  ^  (see  Equation  (9) 

.  *x* 

Chapter  III)  is  equal  to  n  Also  since  B2  R  ,  the  number  of 

inputs  must  equal  the  number  of  outputs  to  satisfy  this  condition.  There¬ 
fore,  in  order  to  maintain  8  inputs,  8  outputs  are  required.  This 
prevented  the  ides  of  measuring  less  outputs  to  place  the  B  matrix  in  the 
proper  format.  Also,  if  n  m  jLt  then  the  C  matrix  is  equal  to  Cj. 

The  first  approach  uses  the  fact  that  C7B7  -  C  B  and  this 
relation  is  used  in  conjunction  with  the  regular  design  procedure. 

Since  C  is  chosen  to  have  full  rank, the  use  of  a  measurement  matrix 
is  not  needed.  However,  B  is  rank  deficient  so  that  KQ  can  not 
be  determined  by  the  relationship 

lo  "  ^2-2)_1  ^  "  <£  I)"'!  “  ®l£l£  (23) 

This  makes  it  necessary  to  use  a  psuedolnverse  for  B  *.  This 

-l 

Inverse  is  noted  as  B^  .  The  psuedolnverse  la  not  guaranteed  to 
have  a  unique  solution. 

Whan  the  B  matrix  psuedolnverse  is  obtained,  the  result  is  a  rank 
•1 

deficient  B^  .  There  are  two  columns  of  saros  (fourth  and* eighth) . 
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Using  the  value  of  Kq  obtained  from  Equation  (23),  and  using 

the  original  C  and  B  matrices,  Equation  (24)  can  be  used  to  see  if 

the  original  £  matrix  can  be  obtained.  The  original  X  matrix  used  to 

find  K  is  X  ■  I.  However,  when  Equation  (24)  is  used,  the 
— o  ~ 

resulting  matrix  obtained  using  TOTAL  (Ref.  8  )  is 


»*o“ 


0.996 

0 

0 

0 

0 

0 

0 

0 


0 

1 

0 

0 

0 

0 

0 

0 


0 

0 

1 

0 

0 

0 

0 

0 


0 

-1 

0 

0 

0 

0 

0 

0 


0.2893 

0 

0.2213 

0 

0.9993 

0.3007 

0.047 

0 


-0.023 

0 

-0.025 

0 

0 

0.9982 

0 

0 


0 

0 

0 

0 

0 

0 

1 

0 


0 

0 

0 

0 

0 

0 

0 

0 


This  is  not  X  that  is  used  to  find  j^.  The  problem  lies  in 
the  psuedoinverse  of  B.  Therefore,  the  psuedoinverse  approach  is 
deemed  unacceptable.  The  next  approach  p>-..s  the  state  and  output 
equations  into  the  form  of  Equations  (9)  and  (10). 

Recombined  Input  Surfaces 

For  the  model  presented  in  Appendix  A,  Equations  (A- 39)  and  (A-40), 
n  ■  8  -  m  -  l .  The  B  matrix  has  two  rows  of  zeros  in  it.  These  two  rows 


of  zeros  cone  from  Che  kinematic  relationships 

9  -  q  (25a) 

•  -  P  (25b) 

If  the  state  equations  are  put  in  the  form  of  Equation  (9),  the  number 
of  outputs  allowed  is  determined  from  the  equation 

n  -  J(  “  number  of  rows  of  zero  (26) 

Since  there  are  two  rows  of  zeros  and  eight  states,  the  number  of 

outputs  needed  is  6.  Since  the  number  of  outputs  must  equal  the  number  of 

inputs,  there  can  only  be  six  Inputs.  Since  there  are  eight  inputs  in  the 

development  of  the  aircraft  equations,  some  of  the  input  surfaces  must  be 

recombined.  The  surfaces  that  are  split  are  the  horizontal  tall  (C  ),  the 

n 

ailerons  ( $ a)  and  the  spoilers  ( % a) . 

The  ailerons  are  combined  so  that  Bj  has  full  rank.  If  this  surface 
were  left  divided,  there  would  be  two  columns  In  the  B  matrix  that  are 
not  Independent.  Combining  the  aileron  input  surface  left  the  model  with 
7  inputs  so  that  one  other  surface  was  recombined. 

The  spoilers  are  chosen  to  be  the  other  surface  to  be  recombined 
for  two  reasons.  First, the  spoilers  normally  act  as  a  unit,  and  second, 
the  horizontal  tall  is  used  in  another  study  as  a  split  surface  (Ref.  16). 

With  only  six  inputs,  the  state  equations  are  rearranged  so  that 
the  B  natrlx  has  the  form  shown  in  Equation  (9) .  The  resulting 


state  equation  matrices  for  Mach  0.6  at  an  altitude  of  15,000  feet  are 


0 

0 

0 

0 

1  0 

0 

0 

0 

0 

0 

0 

0  0 

1 

0 

-32.2 

0 

-0.06006 

-1.897 

0  0 

0 

0 

0 

0 

-0.000676 

-1.006 

1  0 

0 

0 

0 

0 

0.001362 

-4.835 

-34.08  0 

0 

0 

0 

0.05077 

0 

0 

0  -0.01622 

0.05106  -0.8776 

0 

0 

0 

0 

0  -33.43 

-2.443 

0. 7382 

0 

0 

0 

0 

0  5.326 

-8.1 

-32.1 

(27) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-7.55 

-7.55 

0 

-2.78 

12.04 

0 

-0.0675 

-0.0675 

0 

0.02268 

0.2244 

0 

-8.204 

-8.204 

0 

-0.4132 

-1.795 

0.045 A 

0 

0 

-0.01132 

0 

0 

1.053 

0.7866 

-0.7866 

6.744 

0 

0.8808 

5.08 

-0.0385 

0.0385 

0.4656 

0 

-0.0054 

(28) 
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(29) 


The  B2  matrix  has  full  rank  (rank  -^)  and  a  choice  of  outputs 
Is  made  so  that  has  rank  -  Jl  .  The  output  vector  chosen  is 
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This  choice  of  outputs  (where  ~T  is  the  flight  path  angle)  yields  a  C 


matrix 


The  £2  matrix  has  full  rank  (rank  m Jl  )  so  the  first  procedure 
used  with  this  new  model  is  the  unknown  design  procedure 


Unknown  Design 


In  trying  to  use  the  unknown  procedure,  the  primary  criterion  is 


»  —*  '»-r  * 


Chat  rank  G(0)  •  X  •  When  Che  new  model  la  entered  into  the  computer 
program,  the  resulting  G(0)  matrix  la 


G(0)  - 


0 

2.002 

2.002 

0 

-0.3332 

22.48 

0 

5422.0 

5422.0 

0 

436.5 

-1685.0 

0 

0 

0 

0 

0 

0 

0.2781 

0.1524 

-0.1524 

1.98 

0 

0.2203 

-0.4829 

0.01181 

-0.01181 

0.3447 

0 

0.02059 

0 

0 

0 

0 

0 

0 

(33) 


Clearly,  the  rank  of  G(0)  is  not  equal  to  X  •  Therefore,  the 
unknown  plant  approach  cannot  be  used.  The  next  approach  la  to  use  the 
regular  design  since  rank  CjBj  ■  X  • 


Regular  Design 

The  regular  design  produces  a  matrix  according  to  Equation  (13). 

The  necessary  criteria  is  that  both  Cj  and  have  rank  “  J[  .  This 

is  the  case  for  the  outputs  given  in  Equation  (31).  Also,  it  is  necessary 
that  the  transmission  zeros  of  the  plant  be  in  the  unit  disc  (or  in  the 
left  half  s-plane).  The  location  of  the  transmission  seros  is  given  by 
Equation  (14).  Using  this  equation,  it  is  determined  that  the  transmission 
seros  lie  on  the  unit  disc;  not  inside  of  it.  This  corresponds  to  being 
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at  the  origin  in  the  s-plane.  The  result  of  this  is  that  an  uncontroll¬ 
able  node  (  namely  the  forward  perturbation  velocity  u  )  la  present.  The 
reason  the  transmission  zeros  are  on  the  unit  disc  is  that  the  kinematic 
parameters  p  and  q  are  included  in  the  outputs.  Therefore,  to  obtain 
a  stable  response,  these  outputs  cannot  be  controlled.  Removing  them 
makes  Cj  rank  deficient  so  the  necessary  approach  is  to  use  the 
irregular  design. 


Irregular  Design 


A  new  output  vector  is  chosen  so  that  p  and  q  are  not  Included. 
The  resulting  vector  is 


(34) 


Using  this  output  vector  produces  a  C  matrix  that  is 


10  0 


-10  0  0 


0 


c 


0 

1 

0 


0 

0 

0 


1 

0 

0 


0  0  0 


0  0  0  0 
0  0  0  0 
0  0  10 
0  0  0  0 


0 

0 

0 

1 


0 


1  0 


0  0  0  0 


(35) 


The  C  matrix  in  this  case  Is  rank  deficient.  Since  the  B  matrix 

~2  - 

is  not  changed,  has  full  rank.  This  fits  the  form  for  the  ir¬ 

regular  design. 

The  main  objective  of  the  Irregular  design  is  to  choose  a  measure¬ 
ment  matrix  M  that  makes  the  matrix  £2  have  full  rank.  The 
measurement  matrix  M  also  determines  the  location  of  the  transmission 
zeros.  It  is  necessary  to  place  them  within  the  unit  disc  (  or  equiva¬ 
lently  in  the  left  half  s-plane).  The  measurement  matrix  must  have  the 
seme  number  of  rows  as  C^t  and  the  number  of  columns  in  M  is  equal  to 
the  number  of  rows  in  A^.  For  the  case  at  hand,  M  must  have  dimension 
6x2.  The  matrix  £2  can  be  made  to  have  full  rank  by  placing  elements  in  the 
(3,3)  position  and  the  (6,5)  position.  Since  the  matrix  A^  cannot  be 
changed,  the  M  matrix  must  have  elements  in  the  (3,1)  and  (6,2)  locations 
in  order  to  make  C2  full  rank.  The  chosen  M  matrix  is 


✓ 
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(36) 


The  matrices  F^ 

and 

*2 

are 

determined 

by  using  Equation  (16) 

and  (17).  Since 

-  o. 

*1 

is  equal  to 

for 

any  choice  of  M. 

For  the  M  chosen,  F 

2  18 

0 

-l 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

F«  - 

-2 

0 

0 

0.25 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0.25 

0 

The  location  of  the  transmission  zeros  can  be  determined  by  applying 
Equation  (22).  For  a  sampling  period  of  0.01  and  for  the  M  chosen,  the 
two  transmission  zeros  are  located  at  0.96  (or  -3  in  the  s-plane).  This 
location  is  accepted  for  this  design.  A  different  location  might  be 
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better,  end  this  needs  to  be  Investigated  in  future  work. 


A  problem  that  has  been  observed  by  designers  using  high  gain 
techniques  is  that  large  Inputs  are  required  for  "good"  responses.  For 
aircraft  applications,  the  Inputs  have  physical  limits  that  cannot  be 
exceeded.  For  instence,  the  spoilers  cannot  go  negative,  the  flaps  can¬ 
not  go  positive,  and  the  surface  deflections  cannot  exceed  a  set  angle 
limit.  Tracking  simulations  are  used  to  adjust  design  parameters  until 
a  design  is  achieved  that  does  not  exceed  the  aircraft's  physical  limits. 
The  limits  assumed  for  this  study  are  shown  in  Table  II. 


TABLE  II 

Aircraft  Surface  Limits 


SURFACE 

LIMIT 

Rudder 

+ 

30° 

Horisontal  Tell 

+ 

30° 

Aileron 

45° 

Spoiler 

♦ . 

60° 

Flap 

60° 

The  design  is  initiated  with  T  »  0.01,  c<r-  2,  £  -  1,  end  X  »  I. 
It  is  determined  that  ©<r*  1  is  better  than  0<r*  2,  and  that  £  must  be 
0.1.  The  sampling  period  is  satisfactory.  The  main  parameter  that  has  to 
be  adjusted  is  the  I  matrix.  Tracking  Y  produces  inputs  that  are 
too  large  by  a  factor  of  10.  Therefore,  the  first  element  of  the  X 
matrix  is  reduced  by  a  factor  of  10  from  1  to  0.1.  Likewise,  tracking 
fit  r,  and  fl  causes  input  problems  so  the  corresponding  Z  matrix 
elements  have  to  be  reduced.  Figures  3  and  4  show  input  values  obtained 
tracking  fi  before  and  after  adjusting  the  X  matrix  element.  Likewise, 
Figures  5  and  6  show  the  adjustment  needed  for  tracking  r  while  Figures 
7  and  8  show  the  adjustment  needed  when  tracking  fl. 
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Input  (degree*)  I  |  Input  (degree*) 


Figure  3  Input  Response*  Tracking  /&  with  4  th  Element  of  £  Equal  to  1 


Figure  4  Input  Responses  Tracking  with  4th  Element  of  £  Equal  to  0.001 
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Input  (degrees) 


Figure  5  Input  Responses  Tracking  r  with  5  Element  of  l  Equal  to  1 


Figure  6  Input  Responses  Tracking  r  with  5  Element  £  Equal  to  0.01 


Input  (degrees) 


Uuvtfu o  o 


▲  ♦  O  JL  >  u  "  A  W 


. . ~  i. - J. 


Time  (see) 


Figure  7  Input  Responses  Tracking  jj  with  6fc**  Element  of  %  Equal  to  0.01 


Figure  8  Input  Responses  Tracking  ^  with  6cn  Element  of  L  Equal  to  0.001 


After  adjusting  the  design  parameter  so  that  the  inputs  are  not 
bounded, the  final  set  of  design  parameters  are 

T  -  0.01 
C<r-  1 
<£  -  0.1 


I  * 


0 

0 

0 

0.001 

0 

0 


0 

0 

0 

0 

0 

0.001 


The  matrix  is  computed  using  Equation  (21).  The  for  this  design 


-0. 1625E-02 

0. 1767E-03 

-0. 1185E-02 

-0.3826E-02 

-0.5443E-03 

-0.1066E-04 

0.7340E-01 

-0.6928E-02 

0.2528E-01 

0.1062E+00 

0.9722E-02 

0.4446E-03 

-0.7334E-01 

0.9020E-02 

-0.8172E-01 

-0.1062E+00 

-0.9722E-02 

-0.4446E-03 

-0.6519E-02 

0.7085E-03 

-0.4753E-02 

-0.2418E-01 

-0.2183E-02 

-0.4274E-04 

0.3428E+00 

-0.7893E-01 

0.4034E+00 

0. 3999E-14 

0.3661E-15 

0.1674E-16 

-0.7919E-01 

0.8606E-02 

-0.5774E-01 

-0.6333E-15 

-0.5797E-16 

-0.265 IE -17 
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Ay  a  check  on  the  validity  of  thia  K  ,  Equation  (21)  ia  re-arranged 

o 

to  aolve  for  E.  The  reaulting  equation  ia 

"  £  (39) 

The  reault  of  ualng  Equation  (39)  la  the  £  matrix  of  Equation  (38).  Thia 
la  a  necessary  condition,  but  not  a  sufficient  one  to  guarantee  a  good 
design.  The  final  tracker  control  lav  designed  is 

U(kT)  -  Tfcl^  e(kT)  +  T£  *(kT)  (40) 

where 


T  -  0.01 

«5  -  0.1 

C  -  K. 

-o  —1 


The  tracking  responses  and  robustness  of  this  control  law  is  demon¬ 
strated  in  Appendix  B,  and  a  discussion  of  these  results  is  found  in 
Chapter  VI. 


This  chapter  presents  the  different  design  approaches  taken  to 
develop  a  tracker  control  lav  for  the  A-7D  aircraft.  The  design  is 
for  the  flight  condition  of  Mach  0.6  at  an  altitude  of  15,000  feet. 
All  design  procedures  are  used,  but  the  irregular  approach  proves  to  be 
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Che  only  one  that  gives  reasonable  design.  The  sensors  and  actuators 
are  not  used  In  the  simulation,  and  this  is  discussed  in  greater  detail 
in  Chapter  VI.  Responses  obtained  for  the  final  control  law  are  in 
Appendix  B.  All  the  designs  are  accomplished  using  the  computer  package 
discussed  in  Chapter  IV. 


VI.  Results 


Introduction 

This  chapter  discusses  the  results  of  the  tracker  control  law 
designed  in  Chapter  V.  The  simulation  results  obtained  using  this  design 
are  in  Appendix  B.  Each  output  variable  is  tracked  by  Itself  for  a 
total  of  six  different  tracking  commands.  In  each  case,  the  other  five 
outputs  are  conoanded  to  zero.  For  example,  when  the  flight  path  angle 
is  commanded,  the  command  vector  is 


(41) 


In  Equation  (41),  the  tracked  output  is  commanded  to  one  unit,  and 
the  other  outputs  are  commanded  to  zero. 

After  showing  the  tracking  responses  of  the  control  law.  Appendix  B 
presents  the  robustness  results.  First,  the  control  law  is  applied  with¬ 
out  a  rudder  in  the  plant.  Second,  the  control  law  is  applied  at  another 
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flight  condition  (Mach  0.18  at  2,000  feet).  The  robustness  is  discussed 
in  this  chapter. 

To  conclude  this  chapter,  the  effect  of  changing  design  parameters 
(e*r*  T,  £  ,  X)  is  mentioned. 

Tracking  Responses 

The  final  design  presented  in  Chapter  V  does  not  track  lateral  flight 
modes  very  well;  if  at  all.  The  reason  for  this  is  the  fact  that  the  last 
three  elements  of  the  7L  matrix  are  mall.  These  values  had  to  be  made 
small  so  that  inputs  do  not  exceed  their  physical  limits.  The  end  result 
is  that  the  control  law  designed  in  this  thesis  is  a  good  longitudinal 
tracker. 

Looking  at  Figures  B-7,  B-20,  and  B-33  it  can  be  seen  that  "good" 
tracking  is  obtained  for  the  longitudinal  outputs.  What  is  considered 
"good"  is  determined  solely  by  the  designer.  In  all  cases  the  inputs  do 
not  exceed  their  assumed  limits.  However,  some  inputs  appear  to  still  be 
growing  at  the  end  of  the  simulation.  These  inputs  will  stop  increasing 
when  final  steady  state  tracking  is  achieved.  In  one  case  the  spoiler 
is  shown  going  negative  and  the  flap  positive.  This  can  not  be  physi¬ 
cally  done,  but  a  relation  between  the  surface  can  be  derived  that 
gives  the  same  effect  with  opposite  deflections.  Therefore,  this  is 
not  considered  a  problem.  The  interaction  achieved  is  reasonable.  Inter¬ 
action  of  57.  is  obtained  when  tracking  lT  (  see  Figures  B-8  thru  B-12  ). 
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The  Amount  of  interaction  can  be  easily  seen  looking  at  Table  S  where  the 
peak  magnitudes  of  all  the  outputs  are  shown.  By  having  essentially  no 
Interaction,  the  system  responds  as  if  it  is  decoupled.  The  unit  for  the 
inputs  is  degrees.  The  units  for  the  outputs  are: 

1.  For  angles,  the  unit  is  degrees. 

2.  For  rates,  the  unit  is  degrees  per  second. 

3.  For  forward  velocity,  the  unit  is  miles  per  hour. 


TABLE  3 

Figures  of  Merit  for  if  Tracking 


Figure  of  Merit 


Desired  Value 


Achieved  Value 


u  Tracking  Responses 

For  the  case  when  the  forward  perturbation  velocity  is  tracked,  the 
command  vector  is 


v  ■ 


Looking  at  Figures  B-13  thru  B-24,  it  is  seen  that  good  tracking  is  again 
possible.  Table  6  shows  the  desired  figures  of  merits  and  the  achieved 
values  for  the  case  of  u  tracking.  Note  that  the  desired  values  are 
different  for  this  tracking  case  because  faster  response  is  possible. 

Table  7  gives  the  maximum  input  values  needed  for  u  tracking.  Notice 
that  none  of  the  Inputs  exceed  their  limits. 

For  this  case  the  spoiler  is  shown  going  negative  and  the  flap  positive. 
This  can  not  be  physically  done,  but  a  relationship  between  the  surfaces 
can  be  derived  that  gives  the  same  effect  with  deflections  in  the  opposite 
directions.  Therefore,  this  is  not  considered  a  big  problem.  Table  8 
gives  the  maximum  values  of  the  outputs.  By  looking  at  this  data,  it  is 
seen  that  interaction  of  less  than  17.  is  obtained  when  tracking  u. 
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TABLE  6 


Figures  of  Merit  for  u  Tracking 


Figure  of  Merit 

Desired  Value 

Achieved  Value 

CP 

0.75 

0.5 

cs 

1.5 

1.3 

1.15 

1.07 

P 

TABLE  7 

Peak  Inputs  for  u  Tracking  (in  2  sec) 


Input 

Max  Value 

Assumed  Limit 

0.02 

+  30° 

0  x 

$Hr 

-0.735 

t  30° 

K 

0.956 

t  30° 

s. 

0.075 

t  45° 

%• 

-8.4 

+  60° 

Sf 

0.913 

o 

o 

1 

TABLE  8 


Peak  Output  Values  for  u  Tracking 


Output 

Peak  Value 

r 

2.8  E-4 

u 

1.07 

9 

3.2  E-5 

0.0 

r 

0.0 

* 

0.0 

0  Tracking  Responses 

For  Cracking  the  pitch  angle;  the  command  vector  is 


Looking  at  Figures  B-25  thru  B-36  it  is  seen  that  the  tracking  for  0 
is  acceptable  although  not  as  good  as  for  the  other  longitudinal  inodes. 
Table  9  givej  the  desired  figures  of  merit  for  9  tracking  along  with  the 
achieved  values.  The  figures  of  merit  were  achieved  very  well.  Table  10 
shows  the  maximum  input  values  needed  for  9  tracking.  As  in  the  other 
cases,  none  of  the  assumed  limits  are  exceeded,  but  the  spoiler  and  flap 
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try  to- go  negative  and  positive  respectively.  However,  they  both  cross 
the  zero  point  at  the  same  time  so  that  this  problem  could  be  avoided  In 
the  same  way  it  is  for  u  tracking.  The  peak  values  of  the  outputs  are 
given  in  Table  11.  Extreme  interaction  is  obtained  when  tracking  0.  This 
is  due  to  the  physical  relationship  of  output  variables  Y  and  u  to  9. 
When  pitch  angle  9  is  commanded,  there  must  be  a  slow  down,  and  the  air¬ 
craft  must  have  a  small  flight  path  angle.  The  only  way  to  overcome  this 
is  to  introduce  control  surfaces  like  jet  flaps,  and  possibly  canards. 
Physical  reality  causes  the  interaction;  not  bad  designing. 

When  interpreting  the  meaning  of  the  aileron  input  ( &  ) ,  it  must  be 
:  remembered  that  positive  is  when  the  right  aileron  is  up  and  left 

aileron  is  down.  Just  the  opposite  is  true  for  negative  §  .  This  is  the 
case  when  the  ailerons  are  considered  as  one  input. 


TABLE  9 

Figures  fo  Merit  for  9  Tracking 


Figure  of  Merit 

Desired  Value 

Achieved  Value 

t 

2 

1.9 

P 

cs 

2 

1.9 

1.15 

1.10 
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TABLE  10 

Peak  Inputs  for  0  Tracking  (in  2  aac) 


i 

Input 

Max  Value 

i 

Assumed 

I  Limit 

0.17 

+ 

30° 

S  H 

-7.9 

+ 

o 

o 

r 

\ 

-8.0 

+ 

u> 

o 

o 

Sa 

0.69 

+ 

45w 

5s 

-50 

+ 

60° 

• 

_ 1 

7.7 

60° 

TABLE  11 

Peak  Output  Values  for  6  Tracking 


Output 


Peak  Value 


Lateral  Mode  Tracking 

The  lateral  modes  do  not  track  the  command  values  because  the  last 
three  X  matrix  elements  are  so  small.  However,  these  cannot  be  In¬ 
creased  without  exceeding  the  Input  physical  limits  (see  Chapter  V). 

Table  12  gives  the  peak  values  of  the  outputs  being  tracked.  None 
of  the  lateral  modes  even  come  near  the  desired  value  of  1. 


TABLE  12 

Peak  Lateral  Mode  Tracking  Values 


Output  being  Tracked 

Peak  Value 

i  i 

/ 

7.2  E-3 

r 

3.5  E-3 

0 

1.2  E-2 

No  Rudder  Responses 

The  longitudinal  tracking  is  maintained  when  the  rudder  Is  removed 
as  shown  by  Looking  at  Figures  B-73  thru  B-89.  However,  this  does  not  say 
a  lot  for  the  design  since  the  rudder  has  little  input  to  longitudinal  modes 
anyway.  Figures  B-l  and  B-13  demonstrate  this  fact.  A  better  look  at  the 
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capability  of  this  design  to  handle  a  lost  input  as  a  disability  is 
presented  in  Reference  14.  In  Captain  Porter's  thesis,  a  horizontal 
stabilizer  unit  is  removed.  This  has  a  significant  effect  on  the  long¬ 
itudinal  modes.  Table  13  gives  the  achieved  figures  of  merit  without  a 
rudder.  Table  14  gives  the  maximum  input  values  obtained,  and  Table  IS 
shows  the  interaction  when  the  rudder  is  removed.  All  the  results  were 
very  good  with  the  rudder  disabled.  The  only  real  noticeable  difference 
is  that  the  lateral  modes  are  not  zeroed  out  as  well  without  the  rudder 
when  tracking  u.  However,  interaction  is  maintained  at  less  than  17.  for 
this  case. 


TABLE  13 

Figure  of  Merit  without  a  Rudder 


Figure  of  Merit 


Achieved  Value 
Tracking  Y 


Achieved  Value 
Tracking  u 


TABLE  14 


TABLE  15 

Peak  Output  Values  Obtained  without  a  Rudder  Input 


Output 

Peak  Value  Tracking  X 

Peak  Value  Tracking  u 

r 

1.0 

2.8  E-4 

u 

3.4  E-2 

1.07 

9 

4.25  E-2 

3.2  E-5 

0 

0.0 

-4.3  E-4 

r 

0.0 

2.0  E-3 

0 

0.0 

1.5  E-3 
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New  Flight  Condition 

The  tracker  control  law  designed  at  the  flight  condition  of  Mach  0.6 
and  15,000  feet  Is  applied  at  a  different  flight  condition;  l.e.  Mach  0.18 
and  2,000  feet.  In  this  way  the  control  law  Is  checked  for  robustness 
when  faced  with  changing  parameters.  The  control  law  did  not  respond  well 
at  all.  In  fact,  all  of  the  inputs  exceeded  their  limits  after  about  1.5 
seconds,  and  the  commanded  output  grew  unstable  also.  This  trend  is  seen 
in  Figures  B-90  thru  B-96.  This  does  not  mean  that  the  control  theory  does 
not  work,  but  rather  that  this  is  either  a  case  the  theory  does  not  handle 
or  this  Is  a  poor  design.  A  more  robust  design  can  probably  result  if  the 
design  is  done  at  the  flight  condition  of  Mach  0.18.  The  theory  has  been 
shown  effective  on  several  other  plants. 

Parameter  Effects 

In  doing  the  design,  the  effect  of  changing  design  parameters  (cK^,  T, 
£ ,  and  ^  )  is  noticed.  T  acts  like  a  gain  parameter  since  the  control 
law  is  proportioned  to  1/T  or  frequency.  Fast  sampling  corresponds  to  a 
small  T  and  this  acts  like  high  gain.  So  the  result  is  that  decreasing 
T  speeds  up  the  response.  Likewise,  Increasing  T  slows  down  the  system 
performance.  The  parameter  0^r  is  merely  a  direct  change  in  the  damping 
ratio  of  the  system.  Increasing  o<r  produces  less  damping,  and  decreasing 
C*r  gives  more.  The  parameter  £  is  a  normalizing  factor  used  to  scale 
the  outputs  down  to  the  commanded  value.  For  Instance,  with  £  ■  1,  the 
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outputs  for  this  thesis  are  on  the  order  of  10.  Therefore,  £  is  set 
equal  to  0.1  so  the  outputB  can  track  to  a  value  of  1  and  not  10.  The 
parameter  T.  has  diagonal  elements  that  can  be  changed  to  change  the 
magnitude  of  inputs.  Supposedly,  if  the  fifth  input  exceeds  a  desired 
magnitude,  then  the  fifth  £  element  can  be  adjusted  to  reduce  the 
magnitude.  However,  it  is  found  in  this  thesis  that  if  problems  arise 
with  input  magnitude  with  the  fifth  output;  then  the  fifth  element  of  the 
£  matrix  is  adjusted  to  compensate.  This  is  shown  clearly  in  Chapter  V. 
This  is  the  reason  that  this  matrix  is  referred  to  as  an  output  weighting 
matrix  (Ref.  12). 

Summary 

The  sensor  and  actuators  are  not  used  in  the  simulations  because 
acceptable  responses  can  bot  be  achieved.  The  units  are  instead  taken 
out  of  the  design  process  by  approximating  their  transfer  functions  as 
unity.  This  may  be  a  reason  that  the  program  had  difficulty  with  simu¬ 
lations  of  more  than  four  seconds.  A  good  study  on  how  to  include  the 
sensor  and  actuator  models  is  found  in  Joseph  Smyth's  thesis  (Ref.  18). 

This  chapter  discusses  the  tracking  capability  of  the  control  law 
designed  in  this  thesis.  The  result  is  that  longitudinal  tracking  is 
possible,  but  not  lateral.  The  robustness  of  the  design  is  discussed. 

The  design  is  not  robust  when  applied  at  another  flight  condition.  Further 
robustness  checking  of  this  type  is  found  in  Reference  14.  The  ability 
of  the  design  to  handle  a  surface  disability  is  discussed  also.  The  dis- 
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ability  considered  is  the  removal  of  the  rudder.  Since  the  rudder  has 
little  effect  on  longitudinal  motion,  this  disability  is  handled  easily 
by  the  longitudinal  tracker.  The  final  area  presented  in  this  chapter  is 
the  effect  of  design  parameters.  That  section  may  prove  to  be  the  most 
valuable  to  future  designers  using  this  technique.  The  final  chapter 
presents  the  recommendation  of  the  author  for  future  work. 
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VII.  Conclusions  snd  Recommends cions 


Conclusions 

Although  a  control  law  that  can  track  any  of  the  six  outputs  is  not 
designed,  a  fairly  good  longitudinal  tracker  design  is  achieved.  The 
methods  developed  by  Professor  Porter  are  good  in  many  cases,  and  have 
been  proven  so  on  many  examples.  However,  the  irregular  design  technique 
is  actually  a  high  gain  approach  to  the  problem.  It  has  been  noted  by 
several  designers  that  often  high  gain  control  causes  inputs  to  be  too 
large  to  be  useful  in  aircraft  applications.  This  is  the  problem  encountered 
in  this  thesis.  When  the  design  is  adjusted  so  that  inputs  do  not  exceed 
physical  limits,  the  result  is  poor  performance  for  lateral  tracking. 

Perhaps  someone  more  experienced  in  this  type  of  design  can  develop  a 
control  law  that  performs  better.  As  some  engineers  may  describe  it,  this 
requires  more  fine  tuning.  This  thesis  should  provide  a  good  starting 
point  for  a  better  design. 

It  should  also  be  kept  in  mind  that  the  choice  of  outputs  affects  the 
design.  One  choice  of  outputs  places  the  transmission  zeros  on  the  digital 
unit  disc,  and  a  measurement  matrix  can  not  be  used  to  move  them  due  to 
uncontrollability.  A  second  choice  of  outputs  provides  the  capability  to 
choose  a  measurement  matrix  that  can  move  the  transmission  zeros  within  the 
unit  disc.  This  is  the  choice  used  in  this  thesis  design.  A  different  choice 
of  outputs  may  prove  even  better. 
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The  theory  does  not  appear  to  allow  for  the  case  where  the  number  of 
states  equals  the  number  of  inputs.  This  is  the  case  when  all  primary 
control  surfaces  are  split  into  separate  independent  control  surfaces. 
Because  of  this,  it  is  concluded  that  this  is  not  a  good  approach  to  recon 
figurable  control  when  the  plant  is  like  that  found  in  Appendix  A.  Also, 
the  theory  does  not  allow  for  rank  deficiencies  in  the  B  matrix.  This 
prohibits  the  operation  of  the  ailerons  independently. 

The  best  part  ef  the  technique  developed  by  Professor  Porter  is  that 
there  is  no  discretization  process  required.  As  Professor  Porter  puts  it 
the  plant  does  not  change  Just  because  digital  control  is  being  used. 
Digital  concepts  are  inherent  in  the  design  equations,  but  the  plants  used 
are  continuous.  This  makes  the  design  process  much  easier  to  apply  than 
other  techniques. 

Recommendations 

There  are  several  areas  that  need  to  be  investigated  in  the  future. 
The  design  should  be  performed  at  the  flight  condition  of  Mach  0.18  at  an 
altitude  of  2,000  feet.  Work  has  been  done  and  it  appears  that  a  robust 
tracker  can  be  designed  at  this  condition  (Ref.  14).  It  may  be  possible 
to  design  a  more  robust  tracker  around  this  condition.  Also,  a  different 
choice  of  outputs  may  prove  to  be  better.  A  look  at  how  different  output 
choices  affect  the  design  should  be  studied. 

Since  a  longitudinal  tracker  is  designed  in  this  thesis,  future 


research  can  concenCrate  on  designing  a  lateral  tracker.  Then  once  a  good 
lateral  tracker  is  designed,  the  tvo  control  laws  can  be  combined  into  one 
or  stored  separately  in  a  digital  computer.  The  latter  would  probably 
prove  to  be  easier  to  achieve. 

Since  the  computer  program  is  essential  to  design  work,  improvement 
of  MULTI  should  be  performed.  At  present  the  program  is  not  capable  of 
producing  ealcomp  plots.  This  option  needs  to  be  corrected.  Also,  it  is 
desirable  in  many  cases  to  slow  down  a  system  response  by  ramping  or  shaping 
the  inputs  until  they  reach  the  step  value.  MULTI  does  not  have  the  ability 
to  use  ramped  inputs.  Only  step  Inputs  which  are  piecewise  constant  are 
available.  Software  needs  to  be  developed  so  that  the  designer  has  the 
ability  to  use  either  ramped  or  step  inputs. 

The  simulation  package  developed  and  used  by  Professor  Porter  uses  a 
Runge-Kutta  solution  of  the  differential  equations.  MULTI  uses  a  library 
routine  called  01®  (Ref .  11).  It  should  be  determined  in  future  work  if 
using  a  Runge-Kutta  routine  can  give  better  simulation  results  than  ODE, 
and  if  so,  why.  A  final  area  of  needed  improvement  in  MULTI  is  to  increase 
the  total  time  for  a  simulation.  The  simulation  can  not  be  run  for  more 
than  4  seconds  without  running  out  of  CP  time.  Whether,  this  is  due  to  ill 
conditioned  equations,  or  whether  it  is  a  problem  Inherent  in  the  program 
could  not  be  determined.  Research  in  this  area  needs  to  be  done  so  that 
the  simulation  time  can  be  increased. 

While  there  is  much  work  still  to  be  done,  the  most  Important  recom¬ 
mendation  is  that  this  thesis  be  continued  to  gain  a  better  understanding 
of  the  problems  associated  with  high -gain  control  theory  and  the  program 
MULTI. 
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Appendix  A 


Introduction 

This  appendix  presents  three  models  of  the  A-7D  using  linearized 
equations  of  motion.  The  linearized  equations  are  used  to  derive  the 
continuous  state  apace  models  needed  In  the  design  of  the  tracker 
control  law. 

Equations  of  Motion 

The  lateral  and  longitudinal  equations  are  developed  with  coupling 
between  axes  accomplished  by  non-tradltlonal  control  Inputs.  Control 
derivatives  not  normally  associated  with  either  lateral  or  longitudinal 
axes  are  derived  as  described  in  Reference  15.  For  example.  Equation 
(A-l)  is  a  longitudinal  dimensional  control  derivative  produced  by  the 
lateral  rudder  displacement  and  is  proportional  to  the  coupling  coefficient 


r 


SqCp  ^ 


The  linear  equations  assume: 

1.  The  mass  of  the  aircraft  Is  constant. 

2.  The  aircraft  %ody  is  rigid. 

3.  Perturbations  from  equilibrium  are  small. 

4.  The  X,  Y,  and  Z  axes  lie  In  the  plane  of  symmetry,,  and  the 
origin  of  the  axes  Is  at  the  aircraft  center  of  gravity. 

5.  Flow  Is  quasi-steady. 
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6.  The  earth  is  an  Inertial  reference 

7.  Stability  axes  are  used. 

8.  The  aircraft  flies  in  a  straight  and  level  triimned  flight 
condition. 

The  equations  of  notion  are  developed  using  the  inputs  for  the  hori¬ 
zontal  tail  ($H),  the  rudder  (  Sr)»  the  ailerons  (%a),  the  flaps  ($£>, 
and  the  spoilers  (  §  ) .  The  linearized  longitudinal  equations  of  notion  are 
as  follows: 


u  *=  -g  8  cos 


9i  +  x«  “  +  +  X*A  +  \  Sh  *  x  st  S, 

+  xs.Sa  +  x  s.§-  +  x  stsf 


(A-2) 


w  ■  U,q  -  g  9  sin  8-.  +  Zu  u  +  Z_«  +  Z  .  ek  +  Z„  q 

1  °  1  u  ex  of  q 

+  zsh^h  +  zsr^r  +  zOa  +  zO®  +  zsf^f 


(A -3) 


q  -  My  u  +  M .  <X  +  H  q  +  M .  +  M  Sr  +  M  <-  § 

eX  M  «H  ®r  aa 


+  M  +  Ms  Sf  (A-4) 

®  8  °f 


The  lateral  equation  of  notion  are: 

-  -Ujr  +  g  0  cos  9l  +  Y^  /ff  +  Yp  P  +  Yr  r  +  Y%  §H 


+  V*  1  Ys.s-  +  Ys.s-  +  Ysf  Sf 


(A- 6) 


xz 


t  +  +  Lp  p  +  Lr  r  +  LC,t^H  +  LS_^r 


H 


+  L§ft^a  +  L  %%%  +  Lg^ S£ 


(A- 7) 
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IX2 

*  "  —  P  +  V*  +  "p  p  +  Nr  r  +  N  $„  +  N _  5r 

I XX  #H  »r 

+  Ni“aJa  +  *Sa  S‘ “  +  Nff5f  (A‘8> 

0  -  p  +  rtan  »!  (A-9) 

Notice  that  by  using  separate  control  of  the  segments  of  the  lateral 
control  surfaces  (Jr,  Ja  )  it  is  possible  to  excite  longitudinal  motion. 
This  is  developed  later  in  this  appendix.  Also,  note  that  the  longitu¬ 
dinal  control  surfaces  (  Sh  »  S  8  ,  S  f  )  can  excite  lateral  motion. 


Dimensional  Control  Derivative  Equations 

The  equations  for  the  longitudinal  dimensional  control  derivatives 
as  given  in  Reference  16  are  as  follows: 


-  -XC. 


Ds 

x  <Cdu  +  20^) 

Ui 


(A -10) 


(A- 11) 


x«  -  -x  (C^  - 


(A-12) 


where  X  ■  qs/m 


For  the  Z  force  inputs  the  equations  are: 


Zf  -  -zc^  (A-13) 

Z  -  -Z(CLU  +  2CLi) 

Ux  (A- 14) 
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2  -  -Z  (  CL  +  CD,  ) 


where  Z  ■  qs/m 


For  the  moment  Inputs  around  the  y  axis,  the  equations  are: 


"  MC® 

«u  “ 


l  “S 

M  (  Si  +2^m1) 


‘  MCjj, 


M  c  C. 


M.  - 

oC 


M 


m. 


2U. 


M  c 
2U, 


where  M  ■  qsc  /  I. 


yy 


The  equations  for  the  lateral  control  derivatives  are; 


Yr  -  YC. 


Y.  - 


where  Y  “  qa  /  m 


Y  b  Cy_ 


2U, 


Y  b  Cyr 

2U, 
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(A-15) 

(A -16) 

(A-17) 
(A- 18) 

(A -19) 

(A-20) 

(A  -  2 1 ) 
(A-22) 

(A-23) 

(A -24) 


vT 


For  the  L  and  N  moment  inputs  around  the  x  and  z  axes,  respectively, 
the  equations  are  as  follows: 


LP 


L 


r 


N 

N 


§ 

ft 


N 

P 


N 


r 


LC, 

L% 

LbCfp 

2Ui 

LbCir 

2U1 

NC_ 


NC 

n/8 

NbC 


n„ 


2Uj 

NbC, 


nr 


2U, 


(A-25) 

(A-26) 

(A-27) 

(A-28) 
(A -29) 
(A- 30) 

(A-31) 

(A-32) 


where  N  »  qsb  /  lzz  and  L  ■  qsb  /  Ixx 

In  all  of  these  equations,  §  .  ,  $  a  •  $  r  ’  ^s  *  $f  * 

The  data  for  use  in  these  equations  is  given  in  Table  A-l  (Ref.  1  ). 


TABLE  A- 1 


Aircraft  Data  at  Different  Mach  Numbers 


Mach 

0.18 

0.6 

0.8 

Units 

Altitude 

2,000 

15,000 

35,000 

ft 

Weight 

25,338 

25,338 

25,338 

lbs 

Center  of 

gravity 

28.71  7. 

28.71  7. 

28.71  X 

7.  of  m.g.c. 

q 

44.67 

300.88 

435.99 

lbs  /  ft2 

s 

375 

375 

375 

ft2 

b 

38.73 

38.73 

38.73 

ft 

c 

10.84 

10.84 

10.84 

ft 

*xx 

15,475 

15,365 

13,323 

slug- ft2 

X.z 

73,697 

79,005 

79,005 

slug -ft2 

I 

66,566 

69,528 

69,528 

slug -ft2 

yy 

i 

-3870 

-1664 

-2046 

2 

slug 'ft 

xz 
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Derivation  of  Non-Dimensional  Control  Derivative* 

The  control  derivatives  are  found  by  using  the  following  three 
references:  digital  DATCOM  (Ref.  10),  A-7  Aerodynamic  Data  (Ref.  9), 
and  from  Reference  1.  The  DATCOM  program  is  used  to  find  the  non-dimen¬ 
sional  derivatives  which  ere  not  available  from  other  sources.  At  0.8 
Mach,  some  derivatives  are  Interpolated  from  other  flight  condition  data 
with  the  assumption  that  a  linear  relationship  with  Mach  numbers  exists. 

This  is  done  only  in  cases  where  values  can  not  be  derived  exactly.  Table 
A-II  gives  the  derivatives  for  0.18  Mach  at  an  altitude  of  2,000  feet 
assuming  a  landing  configuration.  Table  A-I1X  gives  the  derivatives  for 
0.6  Mach  at  an  altitude  of  15,000  feet  for  a  cruise  configuration.  Table 
A -IV  gives  the  derivatives  for  0.8  Mach  at  an  altitude  of  35,000  feet  for 
a  cruise  configuration.  To  find  the  needed  control  derivatives  using 
digital  DATCOM  (Ref.  10),  the  same  equations  and  methods  developed  by  Potts 
(Ref.  15)  are  employed.  Therefore,  the  reader  la  Instructed  to  see  Appendix 
B  of  Reference  15.  As  a  summary.  Table  A-V  summarizes  all  the  equations 
needed  when  using  DATCOM. 

After  finding  the  needed  non-dimensional  derivatives,  the  dimensional 
control  derivatives  are  found  by  using  the  equations  of  the  previous  section. 
The  dimensional  control  derivatives  for  the  three  flight  conditions  are 
given  in  Tables  A -VI,  A-VIX,  and  A-VXIX. 


units  are  rad"*  ;  only  right  control  surface  values  are  listed  except  for  rudder  and  flaps. 


iraiea  arc  rad"1  ;  only  right  control  surface  values  are  listed  except  for  rudder  and  flaps 


*  traits  are  t*H;  only  right  control  surface  values  are  Hated  except  for  rudder  end  flaps 

#  Indicates  an  interpolated  value 


-47.623  M  0.00049  Y  27.039  Nc  -9.094 


-30.47  M.  0.0057  T  16.206  Mr  0.0 


Continous  Stm  Space  Model 


The  aircraft  model  has  eight  etata  variables  and  is  normally 
considered  to  have  five  control  Inputs.  Theee  are  described  in  Table 
A -IX. 


TABLE  A- DC 

State  Variables  and  Inputs 


u  | 

perturbation  forward  velocity 

cx. 

perturbation  angle  of  attack 

s 

0 

perturbation  aldealip 

T 

A 

q 

perturbation  pitch  rate 

T 

e 

perturbation  pitch  angle 

E 

s 

p 

perturbation  roll  rate 

r 

perturbation  yaw  rate 

perturbation  roll  angle 

horlsontal  stabiliser  deflection 

I 

N 

Sr 

rudder  deflection 

P 

s« 

aileron  deflection 

U 

T 

s. 

spoiler  deflection 

S 

Sf 

flsp  deflection 

The  controls  surfaces  are  split  into  left  and  right  control  surfaces 
(with  the  exception  of  the  rudder  and  flaps) .  In  this  manner  a  model  can  be 


developed  with  8  inputs.  With  each  state  considered  to  be  an  output, 
there  exist  8  outputs  and  thus  a  square  transfer  function  matrix  is  possible. 

In  order  to  derive  a  state  space  model,  eight  equations  involving 
the  time  derivatives  of  u,  c*  ,  q,  6 , /?  »  P»  *.  and  f  are  required.  They 
can  be  obtained  by  using  Equations  (A-2)  through  (A-9)  along  with  the 


assumptions 

that: 

1. 

Xdt*  Z^,  and 

Zq  are  zero 

2. 

U}  -  constant 

and  Vr  -0  -  Wl 

3. 

Stability  axes 

are  used 

4. 

el  -  °  ■  f^l 

5. 

P1  "  Q1  "  R1  “ 

0 

6. 

w  •  Uj«*  and 

*  -  Vf* 

Using  this  procedure  as  outlined  in  Reference  15,  the  equations  for  0, 
o<  »  q*  6,  p  ,  f>,  r,  and  01  can  be  put  in  the  matrix  form 


x 


Ax  +  Bu 


where 


u 

u 

<* 

• 

$r 

q 

q 

e 

e 

x  ■ 

e 

u  • 

S. 

e 

ft 

/ 

e 

S. 

p 

p 

i 

r 

$f 

L jl. 

u*J 

(A-33) 
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0 


g  0 


0 


V 

/Ui 


V»,  1 


0  0 


(Mu  +  >UZ„)  (M  +  MA.)  Mq  0  0 

str  •• 


V 

'  u 

i 


'Vo,  V,  'Vo,  Vo, 


+  MSH 


**lr  +  M^ 

~uT 


M.Zc  +  Ms 

Ot  «■  ia 


M.Z?  +  M 

<*  V, 


y», 


V. 


vh«rt 


i 


# 


Lt  + 


1 


*xr2 

hat  (I**) 


1 


He2 


and  !  represents  fi  t  p,  r,  Sr,  £a,  *nd  $  The  values  of  Li 

/ 

and  Nt  are  given  In  Table  A-X  and  Table  A-XI. 


TABLE  A-X 

B  Matrix  Roll  Coefficients 


0.18  Mach 

m 

0.8  Mach 

,  -7.1425 

,  -33.425 

-69.24 

s 

-1.98 

-2.443 

-1.834 

Lr 

0.8104 

0.7382 

0.6262 

lSh 

-5.1416 

-4.0445 

0.7709 

Hr 

10.556 

5.735 

0.7708 

Hi 

-33.178 

-17.419 

-2.5196 

Ho 

1.2769 

1.858 

0.2675 

/ 

Hf 

7.427 

4.8184 

0.89172  . 

TABLE  A -XI 


B  Matrix  Yaw  Coefficient# 


mam 

0.6  Mach 

0.8  Mach 

V 

1.5614 

5.326 

9.6398 

n'p 

-4.2536 

8.10 

6.492 

Nr 

-18.082 

-32.184 

-39.882 

% 

0.0244 

0.0385 

-0.0034 

NSr 

-9.166 

-5.1952 

-0.8534 

Nsl 

-1.5712 

-0.2328 

-0.0333 

NS» 

0.0908 

0.0997 

0.0376 

0.0284 

-0.0054 

0.00452 

§  »  i  9  I  / 

Note  that  the  values  for  N^,  Njg»  Lj  ,  and  Lg  are  for  the 

right  aide  control  surfaces  only. 

If  the  control  Inputs  are  divided  as  mentioned  previously,  the 


Input  U  vector  becomes 


I 


where  ,  refers  to  Che  right  horizontal  stabiliser  input  and  £ H  refers 

L 

to  the  left  unit.  The  B  matrix  must  also  be  changed  and  the  split  surface 
B  matrix  is  : 


lere  (  ).  indicates  qu 


Using  Che  sects  equations  in  the  fora 


*  -  AX  +  B*U*  (A -36) 

and  substituting  the  values  given  in  this  Appendix, the  desired  aircraft 
aodels  are  constructed  at  the  three  flight  conditions. 

For  the  flight  condition  of  0.18  Mach  at  an  altitude  of  2,000  feet 
(assuaing  a  landing  configuration)  the  A  natrlx  is  as  follows: 
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The  B  aatrix  for  this  condition  la: ■ 


-14.287 

-14.287 

0 

0 

0 

-1.161 

-1.161 

-4.926 

-0.328 

-0.328 

0 

0.2386 

-0.2386 

0.0815 

0.0815 

-0.9573 

-4.774 

-4.774 

0 

1.5685 

-1.5665 

-0.1752 

-0.1752 

-0.1186 

B  - 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.2664 

0.07634 

-0.07634 

0.003281 

-0.003281 

0 

1.0796 

-1.0796 

2.3792 

6.936 

6.936 

-0.2664 

0.2664 

1.558 

-0.0244 

0.0244 

-9.166 

1.5712 

1.5712 

-0.0909 

0.0909 

0.0284 

0 

0 

0 

0 

0 

0 

0 

0 

where  U* 


(A- 38) 


n 


The  alga  convention  eoployed  la:  poaltlve  aad 
edge  down;  poaltlve  ^  f  la  trailing  edge  left;  poaltlve 
up;  poaltlve  %t  la  trailing  edge  up. 


%  *  la  trailing 
% g  !•  epoiler 
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For  tha  flight  condition  of  0.6  Mach  at  an  altitude  of  15 ,000  fact 
(eruiaa  configuration),  the  k  aatrlx  la  aa  follow*: 


The  B  aatrlx  la  as  follows: 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.01622 

0.05106 

-0.8776 

0.05077 

3.426 

-2.443 

0.7382 

0 

5.326 

-8.10  - 

32.184 

0 

0 

1 

0 

0 

94 


r 


t 

I 

-7.55 

-7.55 

0 

0 

0 

-1.39 

-1.39 

-12.04 

i 

i 

-0.0675 

-0.0675 

0 

0.0479 

-0.0479 

0.01134 

0.01134 

0.2244 

C 

-8.204 

-8.204 

0 

2.151 

-2.151 

-0.2066 

-0.2066 

-1.795 

i 

c 

t 

B  - 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.0454 

0.00566 

0.00566 

0.00134 

-0.00134 

0 

* 

c 

0.7866 

-0.7866 

1.0534 

-3.372 

-3.372 

-0.3567 

0.3567 

0.8808 

f 

-0.0385 

0.0385 

-5.08 

-0.2328 

-0.2328 

-0.0997 

0.0997 

-0.0054 

i 

0 

0 

0 

0 

0 

0 

0 

0 

(A -40) 


\ 

t- 

(  The  final  nodel  la  for  a  flight  condition  of  0.8  Mach  at  an  altitude 


of  35,000  feet  (cruise  configuration).  The  A  matrix  for  this  condition 


is  as  follows: 


-0.2872 

-46.418 

0 

-32.2 

0 

0 

0 

0 

-0.0004927 

-1.212 

1 

0 

0 

0 

0 

0 

0.01378 

-2.8197 

-41.176 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.1945 

0.04864 

-0.893 

0.0388 

0 

0 

0 

0 

-69.24 

-1,834 

0.6262 

0 

0 

0 

0 

0 

9.6398 

-6.492 

-39.882 

0 

0 

0 

0 

0 

0 

1 

0 

0 

(4-41) 


c 

95 


Th«  B  nstrlx  for  this  flight  condition  is: 


-1.034 

-1.034 

0 

0 

0 

-0.461 

-0.461 

-0.00816 

-0.00816 

0 

0.00576 

-0.00576 

0.001927 

0.001927 

-0.6034 

-0.6034 

0 

0.15511 

-0.15511 

-0.00624 

-0.00624 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.00537 

0.00074 

-0.00074 

0.000346 

-0.000346 

0.13002 

0.13002 

0.1484 

-0.4239 

0.4239 

-0.0443 

0.0443 

-0.0034 

0.0034 

-0.8534 

-0.0333 

0.0333 

-0.0376 

0.0376 

0 

0 

0 

0 

0 

0 

0 

To  design  s  tracker  control  law,  the  output  equation  has  the  form 

X  m  9* 

For  this  Study  the  desired  comas nd  input  is: 


0 

(A  >42) 


(A -43) 


v  (t)  - 


-  e-o<-ir 

•  <1 

-  e 

-  fi 

-  p 
•  r 

-  0 


(A -44) 


«h«n  T  it  the  flight  path  angle.  The  C  matrix  for  Equation  (A-42)  la: 


In  thia  way  there  are  eight  outpute. 


Senior#  and  Actuators 

A  physical  plant  must  have  some  means  of  sensing  outputs  and  a  means 
to  move,  or  actuate,  control  surfaces.  For  the  A-7D  the  actuators  are 
servos  that  can  be  approximated  by  a  first  order  lag  as  given  in  Reference 
1.  The  transfer  function  of  this  servo  Is: 

I  -  20 

e&  S+20  (A -46) 

where  %  Is  the  deflection  of  the  control  surface  (output),  and  Is  the 
voltage  to  the  servo  (Input);  see  Figure  A-I. 


Figure  A-Z  A-7D  Servo  Block  Diegren 


Zf  Che  Cranefer  function  ia  written  in  atate  apace  form,  the  actuator 
equation  becomes 


%  -  -20  $  +  20  e^ 

Y  -  S 


(A -47) 
(A-48) 


Comparing  this  to  the  standard  state  space  format  it  is  evident  that  for 
the  actuators: 

A  -  -20 
B  -  20 
C  -  1 
D  •  0 


To  measure  the  output  vector  y  ,  three  types  of  sensors  are  assumed. 

An  accelerometer  can  be  used  to  measure  perturbation  forward  velocity  (u). 

The  describing  equation  for  this  case,  found  in  Reference  5  is: 

-  -WjVj  +  Ra  (A -49) 

where  R^  is  the  inertial  acceleration,  is  the  Inertial  velocity 

end  Wj  is  the  rotation  of  the  body  with  respect  to  Inertial  space  (earth) . 

Xf  straight  and  level  flight  is  assumed,  then  •  0  and  Equation  (A-48) 

«• 

becomes  simply: 

*t  -  Rt  (A-50) 
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Sine*  velocity  la  tha  desired  output,  the  output  equation  can  be  written  ea 

y  -  vt  (A- 51) 

froa  Equation  (A-49)  and  (A-50)  it  can  be  aeen  that  for  the  acceleroneter 
that 

A  -  0 
B  -  l 
C  -  1 

0-0 

To  aaaaure  angles,  an  integrating  gyro  can  be  used.  The  transfer 
function  is: 


+  1 


(A-52) 


The  Integrating  gyro  can  be  used  to  Measure  /?  ,  6,  c*r,  0,  and/or  T. 

Typical  values  of  H,  K,  wn,  and  as  given  in  Reference  2  are: 

H  -  10^  gn  •  cs»*/sec 
K  -  3.03  x  10^  ga*cn^/sec 
-  94.25  rad  /sec 
f  -  0.78 

Using  these  values,  the  transfer  function  can  be  written  in  state 
variable  fern  as: 

Ag  +  147.015  A  +  8883.063  Ag  -  293.141  Wt  (A-54) 
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Uttlfli  Ag  ■  Xj  and  Ag 


X2  ,  Equation  (A-52)  can  be  written  as: 


n  .  r  •  ■  1  p'i .  •  1  pq 

[Ijl  -8883.063  -1*7.015  X2  293.1*1  wi 

(A-54) 

Since  Ag  is  the  desired  angle  output,  the  output  equation  can  be  written  as: 


From  Equations  (A-53)  and  (A-54)  it  can  be  seen  that  for  the  integrating 
gyro: 


D  -  0 


Kata  gyro  can  be  used  to  measure  the  rate  output.  The  rate  gyro  can 
measure  q,  p,  and  r.  The  transfer  function  of  the  rate  gyro  is  given 
in  Reference  2  as: 


100 


S  +  1] 


(A-56) 


Type la 1  values  of  H,  Cp,  and  Ip^  as  given  in  Reference  2  are: 


104 

gm 

«  cm2 /sec 

5  x 

103 

gm  •  cm2, 

34 

gm  • 

cm2 

With  these  values,  the  transfer  function  can  be  written  as: 


294.118 

S2  +  147.059  S 


luis  can  be  written  in  phase  variable  form  as: 


i  r°  1  i  r r  ° 

+  W 

x2  0  -147.059  x2  294.188 


(A-57) 


(A-58) 


where 


*,  ■  X, 

ig  -  X, 

Since  Ag  Is  the  desired  output,  the  output  equation  is  the  sane  as 
Equation  (A-54).  prom  Equations  (A-57)  and  (A-54)  it  is  evident  that  for 
the  rate  gyro: 
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Appendix  B 
Simulation  Results 
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Figure  1-9  •  Response  for  #  Tracking 
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Figure  B-12  f  Keeponee  for  r  Trucking 
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Figure  B-18  Sf  Input  for  u  Tracking 
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Figure  B-24  0  Reaponae  for  u  Tracking 
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Figure  B-26  Input  for  9  Tracking 
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Figure  B-27  Input  for  0  Tracking 
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.j.  2  »Jt„— 02-+~ - *• 


Ve  J  Vml^— 

j2-t 

A 

— •  2_>i 

-+■ 

1 

-T 

-+ 

•  I 

“.Ul5 

—r 

I 

-.027 

1 

-1.01 

-+ 

1 

-1.10 

-+ 

I 

-1.37 

1 

-l.ij 

-f 

X 

-1.74 

-1- 

A 


AA 

/'w  * 


L 


.V-  w  . 

V.  vv.  v\ 


.  MW  w'w'rf  w\ 


Time  (eec) 


Figure  B-32  u  Response  for  0  Tracking 
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Figure  B-45  0  Response  for  j  Tracking 
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Figure  B-46  /3  Tracking  Response 
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Figure  B-47  r  Response  for  ft  Tracking 


Figure  B-48  0  Response  for  ft  Tracking 
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Figure  B-53  ^  Input  for  r  Tracking 
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Figure  B-56  u  Response  for  r  Tracking 

136 


I 
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Figure  B-57  0  Response  for  r  Trucking 
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Figure  B-59  r  Tracking  Response 
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Figure  B-60  0  Response  for  r  Tracking 
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Figure  B-62  Input  for  #  Tracking 
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Figure  B-66  %  Input  for  #  Tracking 
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Figure  B-71  r  Response  for  0  Tracking 
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Figure  B-72  0  Tracking  Response 
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hl 


Input  (degrees)  Input  (degrees) 


— .  3  J  » 

«*  3  j>» 

-.b?J 

- .  ->3*i 


— f— 

I 

-+ 

I 

-f 

I 

-+ 


1 

-+ 

1 

-  *- 


— 

1 


k.  W  k 


<V  to  to* to  to  k 


.  V.  wV'm  J\ 

A.  W  *  k  V 

X-^VtoV 


+ 

x 

+ 


+ 

i 


-i- 


TIbb  (see) 


•I - 1  — 

1.  )j 


Figure  B-75  <£ ^  Input  for  f  Tracking  without  a  Rudder 
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Figure  B-90  %  f  Input  for  J''  Tracking  using  Mach  0.6  Control  Law  at  Mach  0.18 
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